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Abstract: We examine the prospects for testing SO(10) Yukawa-unified supersymmetric 
models during the first year of LHC running at ^/s = 7 TeV, assuming integrated luminosity 
values of ~ 0.1-1 fb _1 . We consider two cases: the Higgs splitting (HS) and the D-term 
splitting (DR3) models. Each generically predicts light gluinos and heavy squarks, with 
an inverted scalar mass hierarchy. We hence expect large rates for gluino pair production 
followed by decays to final states with large 6-jet multiplicity. For 0.2 fb _1 of integrated 
luminosity, we find a 5a discovery reach of rrig ~ 400 GeV even if missing transverse energy, 
E™ 1SS , is not a viable cut variable, by examining the multi-6-jet final state. A corroborating 
signal should stand out in the opposite-sign (OS) dimuon channel in the case of the HS 
model; the DR3 model will require higher integrated luminosity to yield a signal in the 
OS dimuon channel. This region may also be probed by the Tevatron with 5-10 fb _1 of 
data, if a corresponding search in the multi-6 + E™ ss channel is performed. With higher 
integrated luminosities of ~ 1 fb , using E™ ss plus a large multiplicity of 6-jets, LHC 
should be able to discover Yukawa-unified SUSY with m g < 630 GeV. Thus, the year 1 
LHC reach for Yukawa-unified SUSY should be enough to either claim a discovery of the 
gluino, or to very nearly rule out this class of models, since higher values of m g lead to 
rather poor Yukawa unification. 
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1. Introduction 



Grand unified theories (GUTs) find a welcome inclusion of supersymmetry (SUSY) into 
their structure in that SUSY tames the gauge hierarchy problem via the well-known can- 
cellation of quadratic divergences [1]. In particular, the GUT group 50(10) is highly 
motivated in that it allows for- in addition to gauge unification- the unification of all the 
matter superfields of each generation into the 16-dimensional spinor representation [2]. The 
matter unification only works if the 15 matter superfields of the Minimal Supersymmet- 
ric Standard Model (MSSM) are augmented by a SM gauge singlet superfield Nf which 
contains a right-hand neutrino (RHN) field. The presence of RHN fields is essential to 
describe data from the past decade on neutrino mass and flavor oscillations; in particular a 
Majorana mass term near the GUT scale, needed to implement see-saw neutrino masses [3], 
should be generated by the breakdown of £0(10) gauge symmetry. In addition to gauge 
and matter unification, in the simplest 50(10) SUSY GUT models- wherein both MSSM 
Higgs doublets reside in a 10 of 5O(10)- one expects Yukawa coupling unification in the 
third generation: f t = f b = f T (= j Vr ) at Mqut- 

Recently, a variety of studies have examined the MSSM(+RHN) to check whether 
the measured values of gauge couplings and third generation fermion masses do indeed 
allow for t — b — t Yukawa coupling unification [4-15]. Essential to the calculation is the 
inclusion of 2-loop renormalization group equations [16] (RGEs) and inclusion of weak scale 
threshold corrections [17] which occur due to the MSSM — ► SM transition in effective field 
theories. These threshold corrections imply that Yukawa coupling unification depends on 
the entire spectrum of SUSY particles, since the SUSY particles enter the various t, b and 
r self-energy diagrams [17]. 

Assuming universal boundary conditions at the GUT scale, the parameter space of 
SO (lO)-motivated SUSY consists of 

mi/2, rn 16 , m 10 , M|>, A , tan/3, sign(fi), (1.1) 

where my 2 is the common gaugino mass at Mgut> m w is the common GUT mass of 
all matter scalars, mio is that of the Higgs soft terms, and Mp parametrizes potential 
splittings in the GUT scale Higgs (and possibly matter scalar) soft terms. Such splittings 
are expected to arise from the breaking of the 5O(10). It has been found that t — b — r 
Yukawa coupling unification can occur in the MSSM within this setup, but only for very 
restricted forms of the soft SUSY breaking (SSB) parameters at Mqut- These include, for 
the case of \i > 0: 

• Al = 2m 2 1Q = 4m? 6 , 

• mie ~ 5 - 15 TeV, 

• mi/2 < mi 6 , 

• tan (3 ~ 50. 
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These boundary conditions were found in Ref. [18] to give rise to an inverted scalar mass 
hierarchy (IMH), wherein first/second generation scalars end up with masses ~ 10 TeV, 
while third generation scalars, Higgs scalars A, H and and \i are of order ~ 1 — 2 TeV. 

A problem with the IMH scheme is that it is inconsistent with radiative electroweak 
symmetry breaking (REWSB), unless the Higgs soft terms are split at Mqut [19]: m "ji u < 
mjj , thus giving rn 2 Hu a head start over tn 2 Hd in its running towards the weak scale. 1 Such 
splitting naturally occurs due to D-term (DT) contributions to all scalar masses arising 
from the breakdown of SO(W). However, applying the splitting to only the Higgs sector 
("just-so" Higgs splitting, HS) 

™% u>d = m 2 w =F 2M 2 D (HS model) (1.2) 

results in better accuracy of Yukawa unification as compared to full DT splitting. Recently, 
it has been shown that DT splitting, combined with the running effect of the neutrino 
Yukawa coupling f Ur and a small mass splitting between first / second versus third generation 
scalars (the DR3 model) can allow for Yukawa coupling unification to a few percent [15]. 
Both the HS and DR3 schemes lead to sparticle mass spectra characterized by 

• m-|(l,2)~10TeV, 

• m~ ^(3) and \i ~ 1 — 3 TeV, 

• m g ~ 300 - 500 GeV, 

• m Wi ^ ~ 100 - 180 GeV, 

• ms ~ 50- 90 GeV. 

Figure [l] shows the location of a large number of Yukawa- unified models in the R vs. m g 
plane, for the HS model (red dots) and the DR3 model (blue dots) obtained through a 
Markov Chain Monte Carlo (MCMC) scan of the parameter space (for details, see [15]). 
Here, the degree of Yukawa unification is quantified as 

max{f t , f b , f T ) n Q \ 

min[ft, fb, It) 

where ft, fb and f T are the top, bottom and tau Yukawa couplings, respectively, evaluated 
at Q = Mqut- In the DR3 case, if we require R < 1.05, then m g < 450 GeV. In the HS 
model, while m g ~ 300 — 500 GeV is favored for low R < 1.05 solutions, it is possible (but 
not likely) to have occassional models with rrig as large as ~ 700 GeV. 

In models with the above listed superpartner spectrum and a bino-like Z\ state, the 
neutralino relic density is computed to be ~ 10 2 — 10 4 times the measured abundance [9,11], 
and the models are seemingly excluded. However, if one invokes the Peccei-Quinn solution 
to the strong CP problem [20-24], then an axion/axino supermultiplet is expected in the 
theory [25]. With an axino of mass nij ~ 1 MeV the neutralinos will decay via Z\ — > aj, 

1 This can be different in non-universal models, see [10, 13, 14]. 
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Figure 1: Scatter plot of Yukawa unified models in the R vs. m g plane, for solutions in the DR3 
model (blue) and the HS model (red). 

which greatly reduces the dark matter density by a large factor: rria/m^. Cold dark 
matter (CDM) solutions can be found consisting of mainly cold axions and thermally 
produced axinos, with a small component of warm axinos arising from Z\ — * ja decay, 
which occurs on time scales of order 1 sec. Since mig ~ 10 TeV, and we expect mig ~ mg, 
the axion/axino CDM scenario allows for a solution to the gravitino BBN problem, and 
can generate re-heat temperatures Tr ~ 10 6 — 10 9 GeV, which can allow for baryogenesis 
mechanisms such as non-thermal [26] or Affleck-Dine [27] leptogenesis to occur [28]. 

Since the value of nig is so low in Yukawa-unified SUSY models, we expect the whole 
scenario to soon be tested at the CERN LHC. 2 LHC has already turned on in Fall, 2009. As 
time progresses, the centre-of-mass energy ^/s will be increased into the ~ 7 TeV regime. 
An integrated luminosity of 0.1 — 1 fb _1 is expected to be collected. 3 Earlier work on 
Yukawa-unified SUSY at LHC with y/s = 14 TeV showed the model to be easily testable 
at LHC [30]. The LHC SUSY events should be characterized by gluino pair production 
followed by three-body decays to states including a high multiplicity of 6-jets. In addition, 
opposite-sign dileptons with mass between 40 — 80 GeV (i.e. between the 7 and Z peaks) 
may be evident. 

In the intervening past year, while LHC recovered from an unfortunate incident involv- 
ing faulty circuits and quenched magnets, the experiments have been measuring millions 
of cosmic muon events. This cosmic data has allowed them to fine-tune their detector re- 
sponse to muons, and to make great strides in alignment of detector elements. We expect 
thus that isolated muons, jets and 6-jets should be readily measurable very early on during 
LHC running, while reliable electron identification (e ID) and even more so reliable E™ 1SS 
measurement, may require additional time to establish. 



2 Indeed experiments at the Fermilab Tevatron collider may also probe up to trig ~ 400 — 430 GeV [29]. 
3 The quoted physics data to be collected at 7 TeV as of November 2009 is 0.2 fb _1 . 
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Figure 2: Total cross-section for gluino pair production with rriq = 10 TeV versus LHC collider 
energy y/s, for m 5 = 300, 400 and 500 GeV. 



In this paper, we expand upon the analysis presented in Ref. [30], and address several 
new issues: 

• We focus on the LHC potential to discover or rule out Yukawa-unified SUSY during 
year 1 of running. 4 To this end, we calculate signal and background production rates 
for the LHC turn-on energy of y/s = 7 TeV rather than the maximal collider energy of 
y/s = 14 TeV used earlier. Part of the effect of LHC turn-on at lower than expected 
energies can be gleaned from Fig. ^ where we plot a(pp — > ggX) vs. collider energy 
y/s, for m-g = 300, 400 and 500 GeV, while taking m- q = 10 TeV. We show both 
LO and NLO QCD results as derived from Prospino [34]. For m g = 400 GeV, LHC 
operating at y/s = 7 TeV yields a cross section of a ~ 10 4 fb. It is expected that, 
after about 0.1 fb -1 of integrated luminosity, the LHC will move up in energy to 
the y/s ~ 10 TeV regime, where a ~ 3 x 10 4 fb. Ultimately, the LHC should move 
up to its design energy of y/s = 14 TeV, where the cross section increases to ~ 10 5 
fb. Various background rates will also change accordingly. In this paper, we take 
a conservative approach, and evaluate all signal and background cross sections at 
y/s = 7 TeV. Increasing the beam energy beyond 7 TeV should only increase the 
SUSY reach projections which we calculate here. 

• We include as well many more background subprocesses than before, including the 
effect of many 2 — ► 3 and 2 — ► 4 body subprocesses. 

• We particularly hone in on what LHC can accomplish with very low integrated lu- 
minosity. After turn-on, some time will be required to examine detector response to 

4 Some additional analyses of early physics prospects at LHC are contained in Refs. [31-33]. 
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well-known SM processes like W, Z and tt production. To be able to use the classic 
SUSY signature of jets + E™ ms production, the measurement of E™ 1SS - which de- 
pends on a knowledge of the entire detector response- will be required. However, in 
Ref. [32,33], it is pointed out that LHC experiments can examine multi-jet + isolated 
multi-muon events in lieu of jets + i?™ lss events as a gain for signal over background. 
We find that for very low integrated luminosity, using either large isolated muon mul- 
tiplicity, or large 6-jet multiplicity, allows m g values of up to 400 GeV to be probed 
with just 0.2 fb _1 of integrated luminosity. (Note that we expect a similar reach for 
the Tevatron in the > 2 - 3 6-jets + Ef iss channel with 5-10 ftr 1 [29].) 

• We also discuss the case when I?™ 1SS measurements and e ID are established. Here we 
find that LHC can explore m g values as high as ~ 630 GeV with 1 fb _1 of integrated 
luminosity. Thus, during year 1 the LHC may well be able to either discover or very 
nearly rule out Yukawa-unified SUSY. 

The paper is organized as follows. We first establish in Sec. ^ two Yukawa-unified 
model lines: one in the HS model and one in the DR3 model. We also examine general 
features in sparticle production and decay for these model lines. In Sec. |3], we present 
some technical details of our signal and background calculations. In Sec. |4], we present 
expectations for early SUSY searches in the multi 6-jets channel 5 without using £™ 1SS cuts. 
We also examine rates for early multi-muon production plus jets without using E™ ss . We 
find that the HS and DR3 models may be distinguishable by measuring the ratio of OS 
dilepton events to multi 6-jet events, since both models produce multi-6-jets at a similar 
rate. However, while OS dimuons from Z2 decay are abundant in the HS model, they are 



relatively scarce in the DR3 model. In Sec. 5.4, we move beyond the 0.1 fb _1 level, and 
calculate the LHC reach for the two model lines using as well E^ BS and e ID for 1 fb _1 
of integrated luminosity. In this case, the 5<r LHC reach should extend to m g ~ 630 GeV, 
enough to cover the bulk of parameter space of these simple Yukawa- unified models. In 
Sec. [7|, we present our conclusions. 

2. HS and DR3 model lines 
2.1 Model lines 



Using the parameter space in Eq. 1.1, Ref. [28] found a large number of SUSY spectral 
solutions with good Yukawa coupling unification in the HS model. We adopt Point B with 
R = 1.02 of this paper as a Yukawa-unified benchmark point, and label it as HSb. The 
HSb input parameters and mass spectra are listed in Table |]. 

To construct a HS model line, we keep most of the above parameters fixed, but allow 
nil 12 to vary. This keeps the Yukawa-unification generally low, but allows us to vary 
nig ~ 3.5m^ ~ 7m^ continuously. We plot the value of R versus m g in Fig. ||. We see 
that at low m g (~ 325 GeV), R < 1.03, while as m g increases, Yukawa unification gets 
worse until m g ~ 700 GeV, where we find R ~ 1.13. 



5 Earlier work emphasizing the utility of the presence of 6-jets in SUSY events was provided in Refs. [30, 
35]. 
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parameter 


HSb 


DR3b 


mi 6 (l,2) 


10000 


11805.6 


mi 6 (3) 


10000 


10840.1 


mio 


12053.5 


13903.3 




nno'7 1 

3287.1 


1850.6 




43.9442 


27.414 


^0 


— 19947.0 


— 22786.2 


tan p 


50.398 


50.002 


R 


1.025 


1.027 




3132. b 


O 1 O O A 

2183.4 


m~ g 


351.2 


321.4 




novo 1 


1 1 Ol A 


tl 


2756.5 


2421.6 


bi 


^77 1 




Tn - 


1 0094 7 


1 1 968 5 


Wi 


116.4 


114.5 




113.8 


114.2 




49.2 


46.5 


m A 


1825.9 


668.3 


m h 


127.8 


128.6 



Table 1: Masses in GeV units and parameters for Yukawa- unified benchmark points HSb [28] and 
DR3b [15]. For the DR3 model, we use M Ns = 10 13 GeV. 

We also adopt from Ref. [15] a DR3 model line, labelled as DR3b, with parameters in 
Table |l[ where mi6(l,2,3), is the scalar mass for the 1st, 2nd and 3th generations and 
Mtv 3 , fu T , A Ut and mjj fl3 are the right-handed neutrino mass, Yukawa coupling, A-term 
and the scalar mass for the sneutrino. We construct a DR3 model line by again keeping 
most parameters fixed, but letting mi/2 to vary. In the DR3 model line, we find R ~ 1.03 
for rrig ~ 325 GeV, while R increases to ~ 1.15 for m g ~ 700 GeV. 6 

Due to the heavy scalar masses in the HS and DR3 model lines, the scalars essentially 
decouple at LHC energies. What results is a low energy effective theory where only g, 
Zi >a and Wf are the new physics matter states. In the next two sections we discuss the 
production cross-sections and decay rates for these states. 

2.2 HS and DR3 production cross sections 

We plot in Fig. |] the leading order gg, W1Z2 and W^W^ production cross sections as a 
function of m g for collider energy ^/s = 7 TeV. From the figure, we see that gluino pair 
production is dominant up to m g ~ 520 GeV (at NLO, it dominates up to m g ~ 560 GeV), 
with cross sections typically greater than 10 3 fb, and in excess of 10 4 fb in the lower gluino 
mass range. Thus, even with integrated luminosities as low as 0.1 fb -1 , we expect hundreds 

6 Although Fig. ^ shows that for some special choice of the parameters lower values of R can be obtained 
for rrig ~ 700, the model lines chosen here represent the general behavior of the bulk of parameter space. 
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of gluino pair events in the upcoming LHC year 1 physics data sample for the HS and DR3 
models. 

Gluino pair production cross sections at the Fermilab Tevatron collider show a large 
increase in rate as rriq increases [29]. This is due to suppression of negative interference 
terms in the qq — * gg subprocess cross section. At the LHC, gluino pair production for 
m g ~ 300 — 500 GeV is dominated instead by the gg — > gg subprocess, which is independent 
of rriq. Thus, gg cross sections show only a slight (~10-20%) increase with increasing rriq 
at the LHC. 

2.3 Sparticle branching fractions in the HS and DR3 model lines 

Since rrig <C rriq in the HS or DR3 model lines, we will get dominant gluino decays into 
three-body modes. The gluino branching ratios will be largely model dependent, but since 
ti and bi are always the lightest squarks, and tan/3 is large [36], the decays will mostly be 
restricted to the following channels: 

• g — > Zi + bb, i = 1, 2 

• g — ► Z\ + tt 

• g -> W^bt or W^bi . 

The general feature m g <C rriq is common to both the HS and DR3 models, since it 
relies mostly on the fact that rri\j2 *C m\§. However, the inclusion of the .D-term splitting 
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Figure 4: Leading order total cross-sections for sparticle production in the HS model (solid) and 
DR3 model (dashed) as a function of the gluino mass for pp collisions at yfs — 7 TeV. The model 
parameters are the same as in Table [j] but with m,xji varying from 30 to 180 GeV. 



for all matter scalars in the DR3 model pushes to lower values, when compared to the 
HS model, where mg £ ~ mi R - 7 As a result we have: 

• DR3: bi ~ 6r and < mg 2 

• HS: b\ ~ bi and ~ mg 2 

Now, since Z 2 is wino-like in both models, it just couples to left-squarks, what suppresses 
the g — > Z2 + bb decay in the DR3 model and favors it in the HS case. This behavior is 
shown in Fig. [|, where the main gluino branching ratios for both models are plotted as 
a function of the gluino mass. From Fig. || it can also be seen that- in the HS model- 
once nig 3> rrit + the g — > W^~bt + c.c. channel starts to dominate, since < 
(for m g > 500 GeV). We can also see that g — > Z\ti becomes relevant for heavy gluinos 
(m g > 600 GeV) and it is enhanced in the DR3 model, where t\ is usually lighter than in 
the HS model. 

From Fig. || we see that, for the HSb benchmark case: 



• BR(g - 


+ Z 2 + bb) = 


63% 


• BR(g - 


^Zi + bb) = 


15% 


• BR(g - 


^Wi + bt) = 


9% 



On the other hand, the DR3b point has: 

7 The stop masses and mixing are basically the same in both models, since the D-term splitting is equal 
for both in and II- 
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Figure 5: Gluino branching ratios for the HS and DR3 model- lines as a function of the gluino 
mass. The model parameters are the same as in Table |l], but with m 1( / 2 varying from 30 to 180 
GeV. 
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Figure 6: Z2 branching ratios for the HS and DR3 model- lines as a function of the gluino mass. 
The model parameters are the same as in Table [| but with 

m i/2 varying from 30 to 180 GeV. 



• BR(g -^Z 2 + bb) = 11% 

• BR(g -^Zx + bb) = 86% 

• BR(g -+Wi + bt) = 0.3% . 

The Z 2 and W\ are expected to decay via three body modes: 

• Z 2 - Z x ff 

• Wf - Zxff' , 

where the decays are dominated by the intermediate virtual W* and Z* diagrams. If 
m g ^ 500 GeV, then the two-body modes W\ — > Z{W and Z 2 — > Z\Z will turn on. 

Putting all segments of the cascade decays together, we expect the HSb signal to be rich 
in 6-jets and opposite-sign/same-flavor (OS/SF) isolated dileptons coming from Z 2 — > Zi££, 
with a small rate of SS dileptons coming from g — > Wxqcf followed by W\ — ► iviZx decay. 
For the DR3b point, we expect the signal to be rich in 6-jets with a harder E™ ms spectrum 
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Gross 


number ol 


SM process 


Generator 


section 


events 


QCD: 2, 3 and 4 jets (p^ > 40 GeV) 


Alp Gen 


3.0 x 10 lb 


13M 


it it i r\ ~i l • j 

tt: tt + 0, 1 and 2 jets 


Alp Gen 


1.6 x 10 lb 


5M 


66: oo + 0, 1 and 2 jets 


AlpGen 


8.8 x 10' fb 


91M 


Z + jets: Z/7(— > It, z^J + U, 1, 2 and 3 jets 


AlpGen 


8.8 x 10 lb 


13M 


VV + jets: vr 3 ^— ► ti>) + 0, 1, 2 and 3 jets 


AlpGen 


1.8 x 10 lb 


19M 


Z + tt: Z/7( — > 11, vv) + it + U, 1 ana 2 jets 


AlpGen 


53 lb 


U.oM 


Z + 66: Z/7(^ i/i/) + 66 + 0, 1 and 2 jets 


AlpGen 


2.6 x 10 3 fb 


0.3M 


W + bb: W ± {-^ Iv) + 66 + 0, 1 and 2 jets 


AlpGen 


6.4 x 10 3 fb 


9M 




MadGraph 


0.6 fb 


1M 


«66 


MadGraph 


1.0 x 10 2 fb 


0.2M 


6666 


MadGraph 


1.1 x 10 4 fb 


0.07M 



Table 2: Background processes included in this study, their cross sections and number of generated 
events. 



(when compared to HSb) due to the direct gluino decay to Zi, but with small rates in the 
multilepton channels. 

3. Event Simulation 

In order to study the discovery potential of the LHC at = 7 TeV we used AlpGen [37] 
and MadGraph [38] to generate the background hard scattering events and Pythia [39] 
for the subsequent showering and hadronization. Table |2| lists the 2 — ► n subprocesses 
included in this study where jets = u, d, s, c and g. For all the processes involving multiple 
jets, the MLM matching algorithm [37] was used to avoid double counting. All the above 
processes were generated at LO, but a relatively low renormalization and factorization scale 
(Q = Vi/6) was used to bring the total cross-sections closer to their NLO values (for more 
details see Ref. [33]). The signal events were generated using Isajet 7.79 [40]. 

A toy detector simulation is then employed with calorimeter cell size Arj x A<j) = 
0.05 x 0.05 and —5 < r\ < 5 . The HCAL (hadronic calorimetry) energy resolution is taken 
to be 80%/^ + 3% for \r]\ < 2.6 and FCAL (forward calorimetry) is 100% /VE + 5% for 
| r/\ > 2.6, where the two terms are combined in quadrature. The ECAL (electromagnetic 
calorimetry) energy resolution is assumed to be 3%o/VE + 0.5%. We use the Isajet jet 
finding algorithm (cone type) to group the hadronic final states into jets. The jets and 
isolated lepton definitions are as follow: 

• Jets are required to have R = \J Aif + A<p 2 < 0.4 and Ex(jet) > 25 GeV. 

• Leptons are considered isolated if they have pr{l) > 5 GeV with visible activity 
within a cone of AR < 0.2 of T,Ef lls < 5 GeV. 

Jets are tagged as 6-jets if they contain a B hadron with Et(B) > 15 GeV, rj(B) < 3 
and AR(B,jet) < 0.5. We assume a tagging efficiency of 60% and light quark and gluon 
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jets can be mis-tagged as a 6-jet with a probability 1/150 for Et < 100 GeV, 1/50 for 
E T > 250 GeV, with a linear interpolation for 100 GeV < E T < 250 GeV (see R. Kadala 
et al. in Ref. [35]). 

4. Early searches without E T nss or electron ID 

During the early stages of data taking at the LHC with integrated luminosity of order ~ 0.1 
fb , it is possible that £?p ltls and electron identification will not be reliable observables 
(see Refs. [32,33]). Therefore we separate our analysis into two stages: 

• early searches, where no E™ 1SS cuts are applied and only muons are considered for 
the leptonic channels, and 

• full analysis, using a minimum £™ 1SS cut and including both e's and /x's. 

In both cases we assume that the 6-jets can be reliably tagged due to their displaced ver- 
tices as reconstructed in the micro-vertex detector or via a non-isolated muon tag. For 
initial searches we apply the following set of minimal cuts, labelled CO: 

CO cuts: 

• Jet cuts: n(jets) > 4 with Et{j) > 50 GeV, rj(j) < 3 and for the hardest jet Ex(jl) > 
100 GeV, 

• Lepton cuts: E T {€) > 10 GeV and r](£) < 2, 

• S T > 0.2, 

• n(6) > 1, 

where St is the transverse sphericity and for now i = u only. 
4.1 Multi 6-jet signal 

As discussed in Sec. |2|the points HSb and DR3b are expected to be rich in 6-jets and possibly 
isolated leptons. In Fig. ^ we plot the 6-jet multiplicity for the signal and backgrounds 
(BG) after applying the CO set of cuts. As expected, the BG distribution falls much faster 
than the signal. Signal exceeds BG for n(b) > 4, where the BG is dominated by multi b 
production (bb and bbbb). We also note that the DR3b case has larger rates for 1 < n(b) < 4 
when compared with the HSb point, due to a lighter gluino. However, the HSb benchmark 
gives a larger signal for n(b) > 5, since in this case Z2 — > Z\bb also contributes, and Z2 are 
produced at large rates from gluino cascade decay. 

We note here that a signal rate above an expected SM BG level may not be sufficient 
to claim a discovery, due to both theoretical and experimental uncertainties in the multi-6 
background rate. However, the overall shape of the n(b) distribution should be to some 
extent self-normalizing, as one can fix the BG levels in the nib) = 0, 1, 2 channels, and 
look for a harder n(b) distribution in the signal case. Thus, some of the uncertainty is 
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Figure 7: 6-jet distribution after CO cuts at the LHC, with y/s = 7 TeV. We show the signal levels 
for the HSb (red) and DR3b (green) points along with various SM backgrounds. 



Results after CO-based selection 




a(n(b) > 3) 


a(n{b) > 4) 


a(OS) 


HSb 


899 fb 


176 fb 


99 fb 


DR3b 


1334 fb 


243 fb 


22 fb 


BG 


1911 fb 


70 fb 


11 fb 



Table 3: Cross-sections for the n(b) > 3, 4 and OS channels after the CO cuts for the points HSb, 
DR3b and the background. 

removed when one looks for an excess in ratios such as a(n(b) = 3)/a(n(b) = 1). Table || 
shows the n(b) > 3 and 4 cross-sections for both points and the background. The signal in 
the n{b) > 4 channel is at the 200 fb level and is well above SM background levels. Such a 
signal may be visible with very low integrated luminosity values ~ 0.05 fb -1 ! 

As mentioned above, a mere excess in one or more of the multi 6-jet channels may not 
be sufficient to claim discovery, due to large uncertainties in the normalization of the high 
jet multiplicities BGs, such as bb, tt and Z + 2 jets. With this in mind, we present some 
signal distributions with distinct shapes from the BG ones, which could help corroborate 
a discovery and provide some information on the sparticle masses. As discussed in Sec. ^ 
events with n(b) > 4 usually come from g — > Zibb decays. Therefore the invariant mass 
of the bb pair is expected to have edges at m g — . This is shown in Figs. || and |||, 
where maxfm^^, m b2 i 2 ] is plotted 8 at the parton level (dashed black line). As expected, 
in the DR3b case the nig — mass edge is much more evident than the rrig — edge, 
due to the large g — > Z\ + bb branching fraction, while the opposite happens for the HSb 
point. At the detector level, the main difficulty in obtaining mox[m 6 ,^,7% 5 ] comes from 

8 Here the index i in bi labels b's coming from the same gluino. 
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Figure 8: Maximum invariant bb mass at parton level for the HSb point without cuts 
(black/dashed) and at detector level for the HSb plus background events (red/solid) with n(b) > 
4, AR(b2,b2) < 1 after the CO cuts (see text). The BG distribution (gray) and the statistical error 
bars for 1 fb _1 of integrated luminosity are also shown. 



combining the correct 6-jets into pairs coming from the same gluino. As pointed out in 
Ref. [30], usually the two hardest 6-jets come from different gluinos. Furthermore, in most 
cases the pair coming from the same g has smaller separation angles. Using these two facts, 
we select b\ and 62 as the hardest and second hardest jets and 62 as the remaining jet that 
minimizes A0(&2, ^2)- Also, to avoid event topologies with no small angular separation, we 
apply the cut AR(p2, ^2) < 1, where R = y A^>(?>2> ^2) 2 + ^(62, ^2) 2 - Using this procedure 
and adding the CO set of cuts we obtain the solid curves shown in Figs. |8| and ||, where 
the BG contribution was added to the signal. We also show the statistical error bars for 
1 fb _1 of integrated luminosity. As can be seen, the invariant-mass distributions have the 
expected shape, although the mass edges seem to require higher integrated luminosity to 
become statistically relevant. 

4.2 Dimuon channels 

In Fig. |l^ we show the muon multiplicity for the HSb and DR3b signal points and the 
BG. For n(fj,) = 0, 1, the BG is well above the signal, but for n(/x) > 2, signal starts to 
dominate over the BG. As expected from the discussion in Sec. g, the DR3b has much 
smaller rates to multileptons. Separating the n(/i) = 2 channel into opposite sign (OS) and 
same sign (SS) muons, we see that almost all the signal comes from the OS dimuon case. 
Due to the small g — * W^tb + c.c. branching fraction, the SS signal is almost 2 orders of 
magnitude below the OS one, which makes it irrelevant for luminosities < 1 fb _1 . As seen 
111 Table g the OS cross-section for the HSb point is around 100 fb and has a discovery 
potential similar to the multi 6-jet channel, while the DR3b benchmark will require more 
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Figure 9: Maximum invariant bb mass at parton level for the DR3b point without cuts 
(black/dashed) and at detector level for the DR3b plus background events (green/solid) with n(b) > 
4, Ai?(&2,fr2) < 1 after the CO cuts (see text). The BG distribution (gray) and the statistical error 
bars for 1 fb _1 of integrated luminosity are also shown. 
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Figure 10: Muon distribution after CO cuts at the LHC, with ^fs = 7 TeV. We show the signal 
levels for the HSb (red) and DR3b (green) points along with various SM backgrounds. In the 
n(/Lt) = 2 bin, the left (pluses) and right (crosses) columns show the background (black) and signal 
components for OS (with invariant mass cuts) and SS dimuons. 



integrated luminosity to be seen in the OS dimuon channel. 

As is well known, the invariant mass of OS muons has a mass edge at — mjf 
since most OS muons come from Z2 — > Zifi + fi~ decays. In Fig. 11, we show the m(/^ + /i _ 



-14- 




m_ -m_ 



BG 

HSb + BG 
DR3b + BG 



50 



100 



150 



m(ii p. ) (GeV) 



Figure 11: OS dimuon invariant mass for the HSb (red) and DR3b (green) points plus background 
events after the CO cuts (see text). The BG distribution (gray) and the statistical error bars for 1 
fb -1 of integrated luminosity are also shown. 



distribution for the HSb and DR3b points. As discussed in Sec. ||, the DR3b point has 
small leptonic rates but may still be visible above background. In both cases, the mass 
edges are visible and give the correct — values. We also show the statistical error 
bars for the combined signal plus background for 1 fb _1 of integrated luminosity. From 



Fig. 11, we can see that the HSb point gives a statistically significant edge while the DR3b 
case may require higher integrated luminosities. 

Finally, we point out that in both the n(b) > 4 and OS dimuon channels the invariant 
mass distributions shown in Figs. ^, ^ and [ll] have distinct features from the BG, which 
should corroborate a discovery claim. 

4.3 Early SUSY search: reach results 

To estimate the LHC reach in the multi-6 and multi-^ channels, we plot in Fig. [l2| the signal 
cross section for the HS and DR3 model lines versus nig using cuts CO plus a). n(b) > 3 
and b). n{b) > 4, along with the expected BG rate. We also plot the 5a discovery lines 
for integrated luminosity values 0.1 and 0.2 fb _1 . The significance in as is derived from 
the p-value corresponding to the number of S+B events in a Poisson distribution with a 
mean that equals to the number of background events [41]. For both the HS and DR3 
model lines, the approximate 5cr LHC reach extends to nig ~ 360 GeV for 0.1 fb _1 , and 
~ 400 GeV for 0.2 fb _1 . We remind the reader that this is comparable to what Tevatron 
experiments can achieve using E™ 1SS cuts and > 5 fb _1 of integrated luminosity [29]. 



The reach using cuts CO and requiring an OS dimuon pair is shown in Fig. 13. In 
this case, the reach in the HS model is similar to the multi-6 reach: LHC should explore 
to m g ~ 360 GeV with 0.1 fb _1 , and m g ~ 400 GeV with 0.2 fb -1 . In the DR3 model 
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Figure 12: Early LHC reach for Yukawa-unified SUSY using cuts CO plus n.j, > 3 and nj > 4. 



line, however, there is no reach in the OS dimuon channel at these low values of integrated 
luminosity. In fact, the rate of multi-6-jet events compared to the rate for OS dimuon 
events would be one way to distinguish early-on whether one might be seeing SUSY in the 
HS or the DR3 model case. 

We also point out that despite giving the maximum reach for both models, the n(b) > 3 
channel has a small signal/BG ratio, what makes it more dependent on the knowledge of 
the background. 
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Figure 13: Early LHC reach for Yukawa-unified SUSY using cuts CO plus requiring OS dimuons. 
5. Analysis including E™ 1SS cut and electron ID 

As the experiments accumulate data, knowledge of the detectors and their response to SM 
background will improve. Also, at some point in time, the LHC center-of-mass energy will 
increase beyond 7 TeV into the 10 TeV range. To be conservative, we will continue our 
analysis assuming ^/i = 7 TeV. Moving to higher values of y/s should only increase the 
possibility of discovering new, high mass matter states. 

As detector response becomes better understood, it will be possible to utilize both 
E™ 1SS and electrons in the analysis. The E™ 1SS variable is well known to be a powerful 
discriminator between SUSY and SM events, and is considered to be the "classic" signature 
for SUSY. In Fig. @ we show the E™ ms distributions for the BG as well as the HSb and 
DR3b points after applying the CO cuts. As expected, the DR3b signal has a harder EJp lss 
spectrum than HSb, due to its large g — > Z\bb branching ratio. 

In order to reduce most of the background and still keep considerably large cross- 
sections for the signal, we henceforth: 

• include the E^ iss > 100 GeV cut and 

• include e's into the CO leptonic cuts 

into our analysis. The cuts CO augmented with the E™ ms cut and inclusion of e's, but with 
no n(b) requirement, will be called CI cuts: 

CI cuts: 

• Jet cuts: n(jets) > 4 with Et(j) > 50 GeV, rj(j) < 3 and for the hardest jet Er(jl) > 
100 GeV, 
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Figure 14: E™ lss distribution for the BG (gray), HSb (red) and DR3b (green) points after the CO 
cuts (see text). 



• Lepton cuts: E T (l) > 10 GeV and r](£) < 2, 

• S T > 0.2, 

• EJp iss > 100 GeV, 
where £ = fi, e. 



5.1 Multi fe-jet + E™ iss channel 

The main effect of adding an E™ ms cut to our previous analysis is the drastic reduction 
of background in the multi 6-jet channel. However, the signal will also be significantly 
reduced, and so will be the statistics in the invariant mass distributions. 

The n{b) distribution after CI cuts is shown in Fig. |l^. Now the signal's peak at 
n(b) =1, 2 is visible above the BG, and the hard distribution in n(b) should be a striking 
signature for both the DR3b and HSb models, since the combined signal plus BG distri- 
bution becomes approximately flat for < n(b) < 2. The S™ 1SS cut reduces most of the bb 
and bbbb backgrounds, leaving tt as the dominant one for n(b) > 1. This results in a consid- 
erable reduction of the BG in the n(b) = 1,2 and 3 bins, where now the signal/background 
ratio is larger than one. The cross-sections for n(b) > 3, 4 are shown in Table |j. Due to the 
large signal cross-section in these bins, the signal could be visible with less than 0.05 fb _1 
of integrated luminosity (in the happy case where an E™ 1SS measurement is immediately 
viable) ! 

The invariant mass distributions from Figs. || and[9| are now re-plotted in Figs. 16 and 



17 after the E™ 1SS cut is included. Now, due to the drastic reduction in the background, the 



m h i distribution is nearly free of BG events. However, the decrease in the signal statistics 
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Figure 15: 6-jet distribution after CI cuts at the LHC, with y/s = 7 TeV. We show the signal 
levels for the HSb (red) and DR3b (green) points along with various SM backgrounds. 



Results after Cl-based selection 




a{n(b) > 3) 


a(n(b) > 4) 


a(OS) 


HSb 


364 fb 


68 fb 


81 fb 


DR3b 


782 fb 


139 fb 


23 fb 


BG 


16 fb 


2 fb 


9 fb 



Table 4: Cross-sections for the n(b) > 3, 4 and OS channels after the CI cuts for the points HSb, 
DR3b and the background. 



makes it impossible to obtain any information on the sparticle masses from the shape of 
max[m bi - bi ,m b2 - b2 }. 

A quantity that may still reveal some information on the gluino mass scale is the 
effective mass, M e s = Yl,PT{jets) + EJp lss , plotted in Fig. |^ for the HSb and DRb points, 
as well as for two points on the HS and DR3 model lines with heavier gluinos (nig = 576 GeV 
and nig = 581 GeV, respectively). 

5.2 Multi-lepton channels 

In Fig. [l9|, we show the isolated lepton multiplicity distribution after cuts CI. In the 
n(£) = channel, the DR3b signal exceeds BG, while in the n(£) = 2 channel, the HSb 
signal exceeds BG. In the n(£) = 1 channel, BG from tt production is larger than both 
signal cases. One can pull out a better signal rate in the 1£ channel by imposing in addition 
a M T (£,E^ iss ) > 100 GeV cut, as is well known [42]. 

We divide the n{£) = 2 channel into opposite-sign/same-flavor (OS/SF) events (e + e~ 
or and same-sign events. For the OS/SF channel, we see the E™ ss cut has little 

effect on the signal-to-BG ratio, since the background, mainly tt production, also has large 
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Figure 16: Same as Fig. |, but including the i^ iss > 100 GeV cut. 
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Figure 17: Same as Fig. |, but including the Ef iss > 100 GeV cut. 



i?™ lss . Again, we expect the entire distribution to be self normalizing, since ti is the 
dominant BG, and one can fix the total ti cross section by normalizing to the 1£ channel. 
Then the n(£) distribution should be harder than expected from just SM physics if a SUSY 
signal is present. 

In Fig. we show the OS/SF dilepton invariant mass distribution after cuts CI. 
Now, despite the negligible background, the reduction in the signal makes the mass edges 
less visible. In particular, the dilepton invariant mass still shows the m^ 2 — edge, 
but with a smaller statistical significance as evident from the error bars. Performing a 
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Figure 18: Effective mass scale of events with n(b) > 3 after CI cuts for the BG, points HSb and 
DRb, and two points on the HS and DR3 model lines with heavier gluinos [m- g — 576 GeV and 
nig = 581 GeV, respectively). 
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Figure 19: Lepton multiplicity distribution after CI cuts at the LHC, with y/s — 7 TeV. We show 
the signal levels for the HSb (red) and DR3b (green) points along with various SM backgrounds. 
In the n(£) — 2 bin, the left (pluses) and right (crosses) columns show the background (black) and 
signal components for OS/SF and SS dileptons. 



different-flavor subtraction may reduce the already negligible BG even further. 
5.3 Jets plus Z -> ££ + Ef iss signal 

We also show in Fig. ^ two additional HS and DR3 cases with nig ~ 525 GeV. In this case, 
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Figure 20: OS dilcpton (/i's and e's) invariant mass for the HSb (red) and DR3b (green) points 
plus background events after the CI cuts (see text). The BG distribution (gray) and the statistical 
error bars for 1 fb _1 of integrated luminosity are also shown. 



according to Fig. ||, the value of is high enough that the two-body decay Z2 — > Z\Z 
is now dominant. We then expect a signature of multiple 6-jets plus E™ ms plus a dilepton 
pair which reconstructs to m(£ + £~) ~ Mz [43]. While signal rates can be very high for 
this channel, SM background is low, coming from processes such as tt and Z + ti. 



5.4 LHC reach for Yukawa- unified SUSY using EJp lss and e ID 

Next, we investigate the full LHC reach for Yukawa-unified SUSY along the HS and DR3 
model lines. First, we require the cut set CI, which includes E™ 1SS > 100 GeV, and then 
require n(b) > 2 or nib) > 3. The SM background level and 5a level for 0.2 and 1 fb _1 



are shown in the plots of Fig. 21, along with expected signal rates from the HS and DR3 
model lines. For the n(b) > 2 case, we find an LHC reach for Yukawa-unifed SUSY out to 
nig = 500 (600) GeV for 0.2 (1) fb _1 . The reach is largely independent of whether one is 
in the HS or the DR3 model. 

For n(b) > 3, the SM background is greatly reduced. In this case, we find a reach 
to m g = 540 (630) GeV for 0.2 (1) fb . Again, the reach in the multi-6-jet channel is 
largely independent of the model line, since both give large numbers of 6-jets + E™ 1SS in 
the final state. We have also calculated the reach in the n(b) > 4 channel; here, the result 
is qualitatively similar to that obtained in the n(b) > 2 case. 



Next, in Fig. |22j, we show the LHC reach for Yukawa-unified SUSY using the CI cuts 
plus requiring a pair of OS/SF dileptons. In this case, the LHC reach is model dependent. 
For the HS model, where we obtain a high rate for g — > bbZ<i decays, we find a reach up 
to m g = 400 (500) GeV for 0.2 (1) fb -1 of integrated luminosity. For the DR3 model line, 
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there is no reach for 0.2 fb , since here the g — ► bbZ\ decay is dominant. With 1 fb , 
however, a signal with 5a significance should be visible for m g ~ 300 — 450 GeV. 

For nig > 450 GeV, the two body decay mode Z2 — > Z\Z opens up and we expect 
to reconstruct Z — ► e + e~ and Z — ► n + u~ within the class of signal events. Here, we will 
adopt cuts CI, but in addition require a OS/SF dilepton pair with 75 GeV< m(£ + £~) < 
105 GeV (cuts CI'). For this topology- > 4 jets +E^ iss + Z -> the dominant SM 

BG comes from ti production. The LHC reach is shown in Fig. 23. As can be seen, no 
significant excess is expected with 0.2 fb -1 of data. However, for > 1 fb -1 , we find that 
this topology can produce a 5a signal for m g ~ 450 — 530 GeV for both the HS and DR3 
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Figure 22: Early LHC reach for Yukawa-unified SUSY using cuts CI for events containing OS/SF 
dileptons. 
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Figure 23: LHC reach for Yukawa- unified SUSY using cuts CV with a reconstructed leptonic Z 
boson. 



model lines. 



5.5 Differentiating the HS and DR3 models 

Here we discuss the possibility of distinguishing the two models discussed so far using 
low luminosities and the channels investigated in the previous sections. As seen from last 
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Figure 24: The ratio of the OS and nib) > 3 cross-sections after CI cuts for the HS and DR3 
model lines plus background as a function of the gluino mass. We also show the statistical error 
bars for 1 fb _1 . 



section results (Figs. 21 and 22) and the discussion in Sec. H for low to moderate m g masses 
we expect the HS and DR3 models to have rather distinct signatures. The first one is rich 
in multi-b jets (n(b) > 4, 5) and OS/SF dileptons coming from the g — > Z2 + bb followed 
by Z2 — > Z/Z* + bb/£ + £~ decays, but has a softer E™ 1SS spectrum due to the two step 
cascade decays. On the other hand, the DR3 model is mainly dominated by g — > Z\ + bb, 
with small cross-sections in the multidepton channels, a moderate number of 6-jets (n(b) > 
3,4) and a harder E™ 1SS spectrum. To explore these features and discuss how well we can 
distinguish both models with year one data, we plot in Fig. the ratio of the OS and 
n(b) > 3 channels. As expected from the above discussion, we see that the ratio is larger 
for the HS model and suppressed in the DR3 case for rrig < 600 GeV. As the gluino mass 
increases, the W\ + bt and Z\ + bb gluino decay channels become available (see Fig. |5|), 
increasing the OS channel in both models. From Fig. ^ we see that the OS/3b ratio should 
be a good discriminator up to m g ~ 450 GeV. For m g > 450 GeV, higher luminosities are 
required in order to distinguish the models. But with high enough luminosities, we should 
be able to tell the models apart up to m g ~ 600 GeV. For even higher m g values, the signal 
features are too similar and the simple OS/3b ratio is no longer useful. 
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6. Comparison with CDF/CMS multijets + E£ iss channel 



Next we would like to compare our results for the multi-6-jets + E™ 1SS signature with a 
multijets + £^? 1SS analysis for SUSY searches proposed and used by CDF [44] and further 
developed in CMS [45]. We implement the following selection. 

• n(jets) > 3 (we take jet pt > 50 GeV), 

• Ef iss > 150 GeV, 

• A(f>(E™ iss ,H T ) < 1, 

• 0.5 < R 2 = jAfiHjl, Ef ss ) + (vr - A0(j2, E™ iss )) 2 < 4, 

• 1.5 < R al = y/ Ef ss )) 2 + (A^(j3, Ef iss )) 2 < 4.25, 

• A<t>(j2,Ef ss ) > 0.35, 
. A(p(ji,E^ iss ) > 0.3, 

• < 1.7 , 

• Finally, we require the highest Et object in each signal event to be hadronic, rather 
than leptonic 9 . 

The angular cuts are introduced in multijets + E™ ms analyses where no lepton information 
is explicitly used in order to guarantee the discrimination of QCD backgrounds which may 
have large E™ ss arising primarily due to jet mismeasurements. In such events, E™ ss typi- 
cally aligns with the 2nd hardest jet, as the hardest jet has a tendency to be mismeasured 
most and eventually becomes the 2nd hardest jet. 

The reach results for CDF/CMS multijets + E™ 1SS cuts with no n(b) requirement can 
be seen in Fig. ^5[ The 5a reach for 0.2 (1) fb _1 of integrated luminosity is found to be 
m-g ~ 370 (430) GeV. 

In Figure |26|, we plot the n(b) signal and background distributions after CDF/CMS 
cuts, which share similar characteristics with the distributions after CI cuts. The dom- 
inance of signal starting with n(b) = 1 illustrates the importance of considering b— jet 
tagging in the multijets + E™ ss analyses. The original CMS analysis made use of the 
B— triggers to significantly enhance the signal/BG ratio, while here we directly cut on 
n(b). We show in Table |5| the cross sections found after implementing CDF/CMS cuts 
with no requirement on n(b) and with n(b) = 1,2,3. 

9 This is an approximation for the leading track isolation step of the indirect lepton veto (ILV) technique 
included in the CDF/CMS selection. Normally in multijets E^ lss analyses where lepton information is 
not explicitly used, W/Z/tt + n jets backgrounds with leptonic W,Z decays can be eliminated by vetoing 
the events whose leading track is isolated. We omit the second step of ILV featuring jet electromagnetic 
and charged fractions since it is designed to eliminate machine and cosmic backgrounds which we do not 
consider here. 
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Figure 25: LHC reach for Yukawa-unified SUSY using CDF/CMS multijets + E^ iss cuts with no 
n(b) requirement. 
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Figure 26: 6-jet distribution after CDF/CMS multijets + E^ iss cuts at the LHC, with y/s = 7 
TeV. We show the signal levels for the HSb (red) and DR3b (green) points along with various SM 
backgrounds. 



As before, we can do much better by requiring a high multiplicity of 6-jets. In Fig. 27 



we adopt the same cuts as in Fig. |25|, but in addition require n(b) > 3. In this case, the BG 
drops by a factor of ~ 150, while the signal drops merely by a factor of ~ 7 for rrig ~ 300 
GeV. The LHC reach increases to m~ g ~ 400 (500) GeV for 0.2 (1) ftr 1 of integrated 
luminosity. 



-27- 



Results after CDF/CMS multijets + £f iss -based selection 




<r(no n(b) req.) 


a(n(b) > 1) 


a(n(b) > 2) 


<r(n(6) > 3) 


HSb 


390 fb 


313 fb 


167 fb 


55 fb 


DR3b 


849 fb 


739 fb 


435 fb 


140 fb 


BG 


1132 fb 


366 fb 


101 fb 


7 fb 



Table 5: Cross-sections for the channels with no n(b) requirement, n(b) > 1, 2 and 3 after the 
CDF/CMS multijets + E^ iss cuts for the points HSb, DR3b and the background. 
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Figure 27: LHC reach for Yukawa-unihed SUSY using CDF/CMS multijets + E^ [ss cuts with 
rib > 3. 

7. Conclusions 

In t—b—T Yukawa-unified SUSY, we expect a characteristic spectrum of superpartners with 
first/second generation squarks and sleptons around 10 TeV, third generation sparticles, 
heavy Higgs bosons and ix around the few TeV level, and very light gauginos, with m g ~ 
300 — 500 GeV (although here we consider even higher values). Thus, at LHC, we expect to 
see gluino pair production at a high rate, followed by gluino decays to bbZi or tbW^ + c.c. 
SUSY searches should therefore exploit the high multiplicity of 6-jets expected in this 
scenario. 

We investigated two model lines- the HS and DR3 cases. The HS case leads to large 
rates for OS/SF dileptons in the final state, while DR3 case does not. We computed 
numerous 2 — ► 2, 2 — » 3 and 2^4 background processes. We found that with just 0.1— 
0.2 fb -1 of integrated luminosity, the LHC discovery reach with yfs = 7 TeV extends out to 
rrig ~ 400 GeV, even without using E™ 1SS cuts. Crucial use is made of the high multiplicity 
of 6-jets in the final state. In the case of the HS model, a corroborating signal appears in 
the n + n~ + jets+ > 1 6-jet channel. 
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C (fb L ) 


0.05 


0.1 


0.2 


1 




HS 


340 GeV 


371 GeV 


400 GeV 


471 GeV 


CO 


Channel 

INI) Q 


n(b) > 3 
340 GeV 


n(b) > 3 
363 GeV 


n(b) > 3 

on A /~1~~\7 

3b)4 GeV 


n(b) > 4 
4b9 GeV 




Channel 


n(o) > 3 


n(b) > 3 


n(b) > 3 


n(o) > 3 




HS 


436 GeV 


480 GeV 


526 GeV 


630 GeV 


CI 


Channel 

T\T> Q 

JJKo 


ra(o) > 3 
4bU GeV 


n(o) > 3 
oUb GeV 


n(b) > 3 
o4o GeV 


n(b) > 3 
boU GeV 




Channel 


n(b) > 3 


n(b) > 3 


n(b) > 3 


> 3 




HS 




341 GeV 


380 GeV 


474 GeV 


CDF/CMS 


Channel 
DR3 


350 GeV 


n{b) > 1 
382 GeV 


n(b) > 1 
420 GeV 


n(b) > 3 
516 GeV 




Channel 


n(b) > 2 


n(b) > 2 


n(b) > 2 


n(b) > 3 



Table 6: Gluino 5a mass reach for the HS and DR3 model lines for different luminosites. We show 
the values for the early search (CO cuts), the full reach (CI cuts) and the CDF/CMS multijets + 
£"™ ISS cuts (CDF/CMS). We also show the channel which optimizes the reach. 

The LHC reach at very low luminosity and without E™ ss is comparable to the Tevatron 
reach in the multi-6 + E™ 1SS channel with 5-10 fb -1 of data. This may lead to a tight 
competition for the discovery or exclusion of the simplest Yukawa-unified SUSY scenario! 

Moving beyond about 0.2 fb -1 , we expect reliable £"™ 1SS resolution and electron iden- 
tification to become available. We find that the LHC reach, using t/s = 7 TeV and 1 fb -1 
of integrated luminosity, will move into the rrig ~ 600 — 650 GeV range for both the HS 
and DR3 model lines, if we require E™ 1SS > 100 GeV along with n{b) > 3. This reach 
is presumably sufficient to rule out Yukawa unification in the DR3 case. In the HS case, 
somewhat larger values of m g can be allowed, although they seem very improbable. The 
LHC reach for 1 fb -1 using the multi-jets + E™ 1SS signature- but making no requirement 
on n(b)- turns out to be much lower. A summary of our various results is presented in a 
convenient form in Table ^. 

At some point LHC energy will be increased to ~ 10 TeV, and this will allow the 
reach to be extended past the values presented here. Thus, our main conclusion is that 
LHC stands an excellent chance to either discover Yukawa-unifed SUSY during year 1 of 
operation, or exclude almost all its model parameter space! 

Last but not least we note that, if the scenario discussed here is realized, a trilepton 
signal from Drell-Yang W1Z2 production should appear when moving into the several fb -1 
regime, giving direct access to the chargino/neutralino sector. 

Note added: While finalizing this manuscript, we found an Atlas note [46] which also 
investigates using multiple 6-jets to enhance the LHC reach for SUSY in the mSUGRA 
model. They examine integrated luminosity values of 0.1 and 1 fb , but take ^fs = 14 
TeV. Their overall results seem quite consistent with the results given here. 



- 29 - 



Acknowledgments 



We thank Harrison Prosper for sharing with us his p-value and significance code. This 
research was supported in part by the U.S. Department of Energy grant number DE-FG- 
97ER41022, by the Fulbright Program and CAPES (Brazilian Federal Agency for Post- 
Graduate Education). The work of SK is supported in part by the French ANR project 
ToolsDMColl, BLAN07-2- 194882. 



References 



[2] 



For recent reviews, see R. Mohapatra, hep-ph/9911272 (1999) and S. Raby, in Rept. Prog. 
Phys. 67 (2004) 755. 

H. Georgi, in Proceedings of the American Institue of Physics, edited by C. Carlson (1974); 
H. Fritzsch and P. Minkowski, Ann. Phys. 93, 193 (1975); M. Gell-Mann, P. Ramond and R. 
Slansky, Rev. Mod. Phys. 50, 721 (1978). 



3] P. Minkowski, phys. Lett. B 67 (1977) 421 



M. Gell-Mann, P. Ramond and R. Slansky, in 
Supergravity, Proceedings of the Workshop, Stony Brook, NY 1979 (North-Holland, 
Amsterdam); T. Yanagida, KEK Report No. 79-18, 1979; R. Mohapatra and G. Senjanovic, 



Phys. Rev. Lett. 44 (1980) 912 



[4] B. Ananthanarayan, G. Lazarides and Q. Shah, Phys. Rev. D 44 (1991) 1613 and Phys. Lett 



B 300 (1993) 245| ; G. Anderson et al. phys. Rev. D 47 (1993) 3702 and Phys. Rev. D 49 
(1994) 3660| ; V. Barger, M. Berger and P. Ohmami, phys. Rev. D 49 (1994) 4908| ; M. 
Carena, M. Olechowski, S. Pokorski and C. Wagner, Ref. [17]; B. Ananthanarayan, Q. Shah 



and X. Wang, Phys. Rev. D 50 (1994) 5980; R. Rattazzi and U. Sarid, Phys. Rev. D 53 



(1996) 1553j ; T. Blazek, M. Carena, S. Raby and C. Wagner, phys. Rev. D 56 (1997) 691E ; 
T. Blazek and S. Raby, phys. Lett. B 392 (1997) 37l|; T. Blazek and S. Raby, phys. Rev. D 



59 (1999) 095002|; T. Blazek, S. Raby and K. Tobe, phys. Rev. D 60 (1999) 11300l| and 



Phys. Rev. D 62 (2000) 05500l|; see also [5]. 



[5] S. Profumo, \Phys. Rev. D 68 (2003) 015006| ; C. Pallis, 
Gomez, G. Lazarides and C. Pallis 



Nucl. Phys. B 678 (2004) 39S ; M 



Phys. Rev. D 61 (2000) 123512| , |^cZ. Phys. B 638 



2002) 165| and phys. Rev. D 67 (2003) 097701^ U. Chattopadhyay, A. Corsetti and P. Nath, 
Phys. Rev. D 66 (2002) 035003|; M. Gomez, T. Ibrahim, P. Nath and S. Skadhauge, \Phys. 



Rev. D 72 (2005) 09500S 



[6] H. Baer, M. Diaz, J. Ferrandis and X. Tata, phys. Rev. D 61 (2000) 11170l|; H. Baer, M. 
Brhlik, M. Diaz, J. Ferrandis, P. Mercadante, P. Quintana and X. Tata, Phys. Rev. D 63| 
(2001) 015007] . 



7] H. Baer and J. Ferrandis, phys. Rev. Lett. 87 (2001) 211805 



T. Blazek, R. Dermisek and S. Raby, phys. Rev. 
Dermisek and S. Raby, 



Lett. 88 (2002) 111804| ; T. Blazek, R 



Phys. Rev. D 65 (2002) 115004; R. Dermisek, S. Raby, L. Roszkowski 



J. High Energy Phys. 0304 (2003) 037; R. Dermisek, S. Raby, L 



and R. Ruiz de Austri, 
Roszkowski and R. Ruiz dc Austri, \,J. High Energy Phys. 0509 (2005) 029| . 



[9] D. Auto, H. Baer, C. Balazs, A. Belyaev, J. Ferrandis and X. Tata, J. High Energy Phys 



0306 2003 023 



[10] C. Balazs and R. Dermisek, J. High Energy Phys. 0803 (2003) 024 



- 30 - 



[11] H. Baer, S. Kraml, S. Sekmen and H. Summy, \j. High Energy Phys. 0803 (2008) 056 . 
[12] W. Altmannshofer, D. Guadagnoli, S. Raby and D. Straub, \Phys. Lett. B 668 (2008) 385 . 
[13] I. Gogoladze, R. Khalid and Q. Shan, \Phys. Rev. D 79 (2009) 115004 
[14] 
[15] 

[16] S. P. Martin and M. Vaughn, \Phys. Rev. D 50 (1994) 2282 
[17] 



D. Guadagnoli, S. Raby and D. M. Straub, \J. High Energy Phys. 0910 (2009) 05£ 
H. Baer, S. Kraml and S. Sekmen, \j. High Energy Phys. 0909 (2009) 005) . 



R. Hcmpfling, \Phys. Rev. D 49 (1994) 6168|; L. J. Hall, R. Rattazzi and U. Sarid, \Phys. Rev 
D 50 (1994) 7048| ; M. Carena etaL, \Nucl. Phys. B 426 (1994) 269| ; D. Pierce, J. Bagger, K. 
Matchev and R. Zhang, \Nucl. Phys. B 491 (1997) 3| . 



18] J. Feng, C. Kolda and N. Polonsky, Nucl Phys. B 546 (1999) 3] J. Bagger, J. Feng and N 



Polonsky, \Nucl. Phys. B 563 (1999) 3j ; J. Bagger, J. Feng, N. Polonsky and R. Zhang, \Phy 



Lett. B 473 (2000) 264. H. Baer,P. Mercadante and X. Tata, Phys. Lett. B 475 (2000) 289; 



H. Baer, C. Balazs, M. Brhlik, P. Mercadante, X. Tata and Y. Wang, fhys. Rev. D 64 
(2001) 015061 . 



[19] 
[20] 
[21] 
[22] 



H. Murayama, M. Olechowski and S. Pokorski, \Phys. Lett. B 371 (1996) 57 



R. Peccei and H. Quinn, \Phys. Rev. Lett. 38 (1977) 1440| and \Phys. Rev. D 16 (1977) 1791 . 
S. Weinberg, \Phys. Rev. Lett. 40 (1978) 223| ; F. Wilczek, \Phys. Rev. Lett. 40 (1978) 279 . 
J. E. Kim, \Phys. Rev. Lett. 43 (1979) 103fc M. A. Shifman, A. Vainstein and V. I. Zakharov, 



Nucl. Phys. B 166 (1980) 493 



M. Dine, W. Fischler and M. Srednicki, \Phys. Lett. B 104 (1981) 199| ; A. P. Zhitnitskii, Sov. 
J. Nucl. 31 (1980) 260. 

For recent reviews on axion physics, see J. E. Kim and G. Carosi, arXiv:0807.3125 (2008); P. 
Sikivie, |hep-ph/0509198| ; M. Turner, \Phys. Rept. 197 (1990)137 . 



Lett. B 138 (1984) 15C| ; J. E. Kim, \Phys. Lett. B 136 (1984) 378 



[26] 



G. Lazarides and Q. Shafi, Phys. Lett. B 258 (1991) 305| ; K. Kumekawa, T. Moroi and T. 
Yanagida, \Prog. Theor. Phys. 92 (1994) 437 ; T. Asaka, K. Hamaguchi, M. Kawasaki and T. 
Yanagida, Phys. Lett. B 464 (1999) 1% 



27] H. Murayama and T. Yanagida, Phys. Lett. B 322 (1994) 349; M. Dine, L. Randall and S 



Thomas, Nucl. Phys. B 458 (1996) 291 



[28] 
[29] 
[30] 
[31] 



H. Baer, M. Haider, S. Kraml, S. Sekmen and H. Summy, JCAP0902 (2009) 002. 
H. Baer, S. Kraml, S. Sekmen and H. Summy, arXiv:0910.2988 (2009). 



H. Baer, S. Kraml, S. Sekmen and H. Summy, \J. High Energy Phys. 0810 (2008) 079 



F. Gianotti and M. Mangano, iep-ph/0505221 ; D. Green, |hep-ph/0601038| H. Baer, V. 
Barger and G. Shaughnessy, \Phys. Rev. D 78 (2008) 095009); J. Hubisz, J. Lykken, M. 



Picrini and M. Spiropulu, \Phys. Rev. D 78 (2008) 075008| ; M. Mangano, \Eur. Phys. J. C 59 
(2009) 373 . 



[32] H. Baer, H. Prosper and H. Summy, \Phys. Rev. D 77 (2008) 055017] ; H. Baer, A. Lessa and 
H. Summy, \Phys. Lett. B 674 (2009) 49| ; J. Edsjo, E. Lundstrom, S. Rydbeck and J. Sjolin, 
arXiv:0910.1106 (2009). 



- 31 - 



H. Baer, V. Barger, A. Lessa and X. Tata, |j. High Energy Phys. 0909 (2009) 063 
W. Beenakker, R. Hopker, M. Spira, Eiep-ph/9611232j (1996). 



U. Chattoppadhyay, A. Datta, A. Datta, A. Datta and D.P. Roy, \Phys. Lett. B 493 (2000) 



127|; P. Mercadante, K. Mizukohi a nd X. Tata, \Phys. Rev. D 72 (30005) 035009|; S. P. Das et 



al. Eur. Phys. J. C 54 (2008) 645; R. Kadala et al. Eur. Phys. J. C 56 (2008) 511 



H. Baer, C. H. Chen, M. Drees, F. Paige and X. Tata, Phys. Rev. Lett. 79 (1997) 986 and 



Phys. Rev. D 58 (1998) 075008 



M. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. Polosa, J. High Energy Phys. 0307 



(2003) 001 



T. Stelzer and W. F. Long, \Comput. Phys. Commun. 81 (1994) 357 . 



T. Sjostrand, S. Mrenna and P. Skands, J. High Energy Phys. 0605 (2006) 026 



ISAJET v7.79, by H. Baer, F. Paige, S. Protopopescu and X. Tata, hep-ph/0312045; for 



details on the Isajet spectrum calculation, see H. Baer, J. Ferrandis, S. Kraml and W. Porod, 
Phys. Rev. D 73 (2006) 015010| . 



H. B. Prosper, private code for calculation of significance from p-values. 



H. Baer, C. H. Chen, F. Paige and X. Tata, \Phys. Rev. D 53 (1996) 6241 



H. Baer, X. Tata and J. Woodside, \Phys. Rev. D 42 (1990) 1450] H. Baer, C. B alazs, A. 
Belyaev, T. Krupovnickas and X. Tata, \.J. High Energy Phys. 0306 (2003) 054 . 



M. Spiropulu, Ph. D. thesis, UMI-99-88600; A. A. Affolder et al. (CDF Collaboration), Phys 



Rev. Lett. 88 (2002) 04180l||hep-ex/010600l| ; S. Sekmen, Ph. D. thesis, CMS TS-2009/025 
S. Sekmen, Ph. D. thesis, CMS TS-2009/025 

Discovery potential for super symmetry with b-jet final states with the Atlas detector, (The 
Atlas Collaboration), Atlas note ATL-PHYS-PUB-2009-075 (2009). 



- 32 - 



